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We analyze single-particle electronic and two-particle magnetic properties of the Hubbard model 
in the underdoped and optimally-doped regime of YBa2Cu307_5by means of a modified version of 
. the fluctuation-exchange approximation, which only includes particle-hole fluctuations. Comparison 

of our results with Quantum-Monte Carlo (QMC) calculations at relatively high temperatures (T ~ 
1000-R") suggests to introduce a temperature renormalization in order to improve the agreement 
between the two methods at intermediate and large values of the interaction U. We evaluate the 

Q temperature dependence of the spin-lattice relaxation time Ti and of the spin-echo decay time 
T2G and compare it with the results of NMR measurements on an underdoped and an optimally 
£ — ■ doped YBa2Cu307_isample. For U/t — 4.5 it is possible to consistently adjust the parameters 

of the Hubbard model in order to have a good semi-quantitative description of this temperature 
dependence for temperatures larger than the spin gap as obtained from NMR measurements. We 
D ' also discuss the case U/t ~ 8, which is more appropriate to describe magnetic and single-particle 

5^ , properties close to half-filling. However, for this larger value of U/t the agreement with QMC as 

■ well as with experiments at finite doping is less satisfactory. 

^ ; PACS numbers: 71.27.+a, 74.72.-h 76.60.-k 
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I. INTRODUCTION 



Intensive research of the last decade made clear that antifcrromagnctism and superconductivity are the two domi- 
nating properties of high-temperature superconductors. Indeed, the fact that these two states of matter do not exclude 
each other and that their fluctuations coexist in an extended parameter range suggests a close relation between them. 
This has been the main motivation for the recently proposed S0(5) theory of superconductivity which unifies antifer- 
romagnetism and superconductivity on the basis of a common symmetry principle. Em Here, as well as in more phe- 
nomenological approaches to the high-T c compoundsuu which relate their underdoped properties to remnants of the 
antiferromagnetic order, the key to understanding the driving mechanism behind superconducting pairing lies in their 
■ magnetic properties. The minimal microscopic model considered to describe strong correlations effects is the Hubbard 
model with nearest-neighbor hopping t and pcwsjte repulsion U. Using Quantum-Monte-Carlo (QMC) simulations 
ON ■ combined with Maximum Entropy techniquesEJO this model has recently been shown to reproduce salient features m 
the underdoped photoemission experiment in particular the pseudogap and its doping and momentum dependence.t3 
The latter has uniquely been related to the momentum and doping dependence of magnetic correlations. c3 In order to 
describe more accuratebz-.the Fermi surface of the cuprate materials, which appears to be closed around the antifer- 
romagnetic poiat (tt, n)f3 as shown in several high-resolution angular-resaUsed-photoemission spectroscopy (ARRES) 
measurement^ on Sr 2 Cu0 2 Cl 2 H, YBa 2 Cu 3 07-5 (YBCO) with S « O.lB and Bi 2 Sr 2 CaCu 2 8+l5 (Bi2212),BB it 
has been suggested to extend the model by introducing additional longer-range hoppings, namely second (t 1 ) and 
O ■ third {t") neareal n£ighbors. This has been enforced by comparison of ARPES data with numerical calculations on 
the t - J modelErall3. Qualitatively, a large range of values of t' and t" yield a Fermi-surface with the appropriate 
shape. However, it is pertinent to specify more precisely the values of the parameters t' and t" which give simulta- 
neously a good qualitative description of other properties, in particular magnetism. In the present work, we combine 
two different many-body techniques, namely, QMC, and a modified version of the FLEX approximation^ (here re- 
ferred to asJVIFLEX, for clarity) whereby particle-particle fluctuations, as well as vertex corrections in two-particle 
correlationsc3 are neglected. The two techniques have advantages and disadvantages in different parameter regimes 
(Coulomb correlation U/t, temperature T, system size N). Here, we intend to provide a definite link between single- 
particle, i.e. photoemission (ARPES) and two-particle, i.e. magnetic excitations. This link may be a useful guide and 
serve as an input not only in the unifying SO(iLtheory but also in phenomenological constructs such as the nearly 
antiferromagnetic Fermi liquid theory (NAFL).trul In practice, we carry out a diagrammatic, self-consistent study of 
single- and two-particle response functions of the Hubbard model with inclusion of longer-range hopping terms up 
to third neighbors t' and t". Our aim is to find a reasonable set of parameters for the model, which consistently 
describe at the same time magnetic (NMR) and electronic (ARPES) properties of doped YBa 2 Cu307_«5 compounds, 
in particular, the Fermi-surface, the band dispersion, the spin-lattice relaxation time T% and the spin-echo decay time 
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T2G- More specifically, we want to adjust these parameters, in particular t' and t", such that the magnetic properties 
are reproduced at least in a semi-quantitative way (and not just qualitatively) at least within a reasonable error. We 
will show indeed that a careful tuning of the parameters is important since a change of t' and t" by only 25% changes 
the result for the T± and Tiq by 100% or more. Since we want to compare theoretical and experimental results at 
moderate antiferromagnetic correlation lengths, it is important to perform the numerical calculations at large system 
sizes at such low temperatures, which at present are not accessible by Quantum Monte-Carlo calculations (at least for 
dynamical correlation functions). For this reason, we will use a refined diagrammatic technique (MFLEX), whereby 
particle- hole diagrams with self-consistently determined Greens functions are summed (See, e.g., Ref. |2?]), which al- 
lows us to work on 64 x 64 (and even larger) lattices down to temperatures T ~ t/50. Since we want to carry out the 
calculation with values of the interaction of the order of the system's bandwith, where perturbational approaches are 
uncontrolled, it is important to compare our results with Quantum-Monte-Carlo (QMC) calculations, which provide 
essentially exact results, in the temperature range accessible to this method. We will show that up to intermediate 
values of U (U/t < 6) our diagrammatic results agree quite well with QMC provided one allows for a renormalization 
of the temperature T. Our idea thus amounts to use the MFLEX calculation to extrapolate QMC calculations to low 
temperatures and large system sizes, which are not reachable by QMC simulations. 

The values of t' and t" giving the best agreement with NMR results turn out to depend on U/t. Consistent results 
are obtained for U/t = 4.5, t = 250meV, t'/t = -0.2, and t" /t = 0.15. For U/t = 8 we need a larger \t' /t\, namely 
t'/t = —0.4, although the comparison with experiments is less satisfactory in this case, possibly due to the fact that 
our diagrammatic calculation is less reliable for large U/t. These values of t' and t", especially the ones obtained 
for U/t = 8 are quite similar tp^thc-Qnes obtained by several authors who compared ARPES data with numerical 
calculations on the t — J modeoua. 

Our paper consists of two main parts. In the first part we compare the diagrammatic results with QMC, while in the 
second part we describe the NMR experiments. More specifically, in Sec. [Il] we shortly introduce the Hubbard model 
and describe the MFLEX approximation. We then discuss single-particle properties like the Fermi-surface and the 
quasiparticle dispersion in Sec. Ill, followed by a detailed comparison of MFLEX and QMC results including single- 
and two-particle properties in Sec. [V. In Sec. [v[ we justify our choice of the hopping parameters by comparing our 
numerical results with NMR data on YBa2Cu3 07_a. Finally, we summarize and draw our conclusions in Sec. N/u. 



II. MODEL AND TECHNIQUE 

The Hamiltonian of the Hubbard model is given by: 

H = X!( £k ~~ V) c l,* c k,* + U^2n t ^n hl , (1) 

k.cr i 

where the bare energy dispersion 

ek = — 2t (cos k x + cos k y ) — At 1 cos k x cos k y 

-It" [cos(2fc x ) + cos(2fc y )] (2) 

includes nearest-neighbor (t) and longer-range hopping processes (t',t"). Here, Ck CT (c k(T ) annihilates (creates) an 
electron with momentum k and spin a. The chemical potential fi adjusts the mean particle number (n) with the 
doping x so that (n) = 1 — x. 

The fluctuation exchange approximation includes the interaction of the electrons with density, spin, and—pairing 
fluctuations in infinite order. According to the conserving approximation scheme in the Baym and Kadanof£3 sense, 
the self-energy is obtained by differentiating an approximate generating functional with respect to the full Green's 
function. For the approximation to be conserving in the two-particle channel it is also necessary to calculate the 
two-particle interaction by taking the second functional derivative of the same generating functional with respect to 
the Green.'s function. This leads to a rather complicated set of coupled equations which can be solved only on small 
systemsBj Within the MFLEX-iapproximation, the electron self-energy S evaluated in the imaginary (Matsubara) 
frequency representation readscJ: 

E(k, iw m ) = — ^2 il/ n)G(k - q, iuj m - iv n ) (3) 

where T is the temperature, N the system size and V(q 7 iv n ) the effective interaction resulting from a geometric series 
of bubble and ladder diagrams 
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V(q,iv n ) = C/ 2 Qxs(q,^n) + ^Xc(q,w n ) - xfawSj (4) 

with 

X(q, w«) = --^ 2J G ( k + <1j iw ™ + ^n)G(k, iuj m ) (5) 

Xc/s(q,^n) = 7X77 7 ( 6 

l±Z7x(q,^n) 

With respect to the FLEX approximation, the MFLEX neglects particle-particle fluctuations which turn out to 
be of minor importance in the parameter range we are considering, i. e.,,-Glose to the antiferromagnetic instability 
where the effectiveunteraction is dominated by the spin- fluctuation part This approximation is conserving at the 
one-particle levelEj, since the self-energy is obtained as a functional derivative of a free energy functional, containing 
particle-hole fluctuations only, with respect to the Green's function. However, the procedure to obtain two-particle 
correlation functions in a conserving way is much more complicated numerically, since functional differentiation of 
the self-energy with respect to the Green's-function would include, beyond the standard ladder and RPA diagrams 
we are considering, also vertex correctionscS. In the present work, we will neglect these vertex corrections and thus 
our MFLEX approximation is not conserving at the two-particle level. This allows us to evaluates physical quantities 
at real frequencies for larger system sizes and smaller temperatures allowing comparison with experimental results. 
Moreover, according to Dahm and TewordtEj, the corrections coming from these additional diagrams seem to be small 
in a similar parameter range. Since we are mainly interested in spectral densities of one- and two-particle correlation 
functions at finite frequencies we should-eventually analytically continue Eq. to real frequencies by inverting the 

corresponding Laplace transformation.^ This inversion, however, introduces errors for large real frequencies w > t 
due to the, exponential kernel of the transformation. In order to avoid these uncertainties we will employ a recent 
approachl!3 which deforms the frequency sums (Eq. [|||) to the line u> + iS with 5 < nT/2 close to the real axis and 
carry out the self-consistent calculation directly on this line. From this line, we can continue analytically our results 
to the real axis with a much better accuracy, since the imaginary part of the susceptibilities at to + iS already shows 
most of the features that are present in the true spectral functions on the real axis. 

The magnetic properties of the Hubbard model are related in linear response to the retarded spin-spin correlation 
function 

/•oo 

- Xzz {q,co) = -i dt e i " t ([S*(q,t)S*(-q ) 0)]> (7) 
Jo 

with S z (q) = l/2Nj2 e ~ lIli ' i ( n i,1 ~ n i.i)- Fox_simplicity, and in order to achieve larger system sizes and lower 
temperatures, we will neglect vertex correctionso in the coupled MFLEX equation for two-particle correlations thus 
calculating the spin response function with the "bubble" sum with the dressed Green's functions obtained within the 
MFLEX formalism and use 

X«(q,w) = - — -. 8 

1 - c^x(q,^) 

Here, x(q, u>) is given by Eq. (||) after continuation to the real-frequency axis. Neglecting these diagrams, makes the 
approximation not conserving, as discussed above. 



III. SINGLE-PARTICLE PROPERTIES 



We focus our study on the properties of two YBa2Cu307_<5 samples, an underdoped one with 6 = 0.37 and a nearly 
optimally doped one with 5 = 0. For the sake of our comparison, however, we need to know the appropriate value 
of the hole doping x to use in our Hamiltonian associated with these two oxygen concentrations. The question of 
how many holes go into the Cu02 layers for a given oxygen content in YBa2Cu30y_5 is quite controversial. Presland 
et alE3 have suggested the empirical formula 

T c /T c , max = l-82.6(z-0.16) 2 (9) 

which relates the doping x to the ratio of T c /T Cimax , where T Cimax is the critical temperature at optimal doping. This 
relation is particular useful in conjunction with measurements of the thermoelectric power at room temperature, since 
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this quantity shows a generic dependence on the hole concentration.E-aH Thus, a measurement of the thermoelectric 
power at room temperature is uniquely related to the doping concentration. According to Eq. (^J), we estimate 
x m 0.08 for the underdoped YBa2Cu306.63 sample and x « 0.16 for the fully oxygenated YBa2Cii307. For these 
doping levels, the parameter set t' = —0.2t and t" — 0.15i yields a Fermi-surface in good qualitative agreement with 
experiments in the sense that it is closed around (jr, tt) and shows a large curvature, see Fig. [I]. Alternative parameter 
sets yielding a similar Fermi-surface are, e. g., t' = — 0.45i,i" = and t' = —0.38t,t" = 0.06i. Of course, for an 
interacting system one expects this bare Fermi-surface to change, when the interaction U is turned on. Indeed, in 
Fig. ^ we show the Fermi-surface for various values of U, including U = 0, obtained by our self-consistent MFLEX 
calculation at the temperature T — 0.02£. Here, the Fermi-surface is defined by the k points matching the condition 

e k + Rc E(k,cj = 0) - (i = . (10) 

Notice that the Fermi-surface is unambiguously defined only for T — 0. At finite T, other alternative definitions (like, 
e. g., the local maximum of the spectral function at lo = 0, or the local maximum of |Vk(^k}|) may give results 
differing an amount of order T jvp from Eq. (|ic|). As one can see, the interaction U modifies the Fermi surface, as 
expected, especially in the regions close to (tt,0) and (tt/2,tt/2) [and symmetrically related points]. By increasing 
the interaction from U/t = to U/t = 8 the curvature of the Fermi-surface is smoothly reduced. The effect of the 
interaction is thus to increase the .magnetic fluctuations by pushing the Fermi-surface closer to nesting with a wave 
vector Q equal or close to (tt, 7r).E3E3 On the other hand, the decrease with increasing U of the area inside the Fermi- 
surface in the region close to (0, tt) seems to be compensated by its increase in the region close to (tt/2, 7r/2). Whether 
this compensation is exact as suggested by the Luttinger theorem is not clear ,t3 since it is difficult to extrapolate the. 
result to T — > were the Fermi-surface would be well defined. The strong effect of U on the quasiparticle dispersior£3 
near the point (tt, 0) is shown in Fig. ^[ The main effects of the interaction are (i) to flatten the dispersion near (tt, 0) 
and (ii) to decrease the binding energy of the quasiparticles Ax = E(tt, 0) — [i close to this pointo. Specifically, the 
energy scale Ax is seen to change from being of the order of the bandwith Ax ~ 0.8t at U = to the order of the 
magnetic excitation J — 4t 2 /U at finite U (more precisely, we have Ax roughly 0.15 • J). The effect of U is thus to 
pin the quasiparticle dispersion at (tt, 0) to the chemical potential. This pinning is related to the onset of a pseudogap 
as was already pointed out in Refs. |l3],[40| Nevertheless, the flat region approaches the Fermi surface with increasing 
hole density. 



IV. COMPARISON OF MFLEX WITH QMC 



Before comparing the magnetic properties obtained in MFLEX and QMC, we first analyze the single-particle 
spectral function A(k,oj) = — ImG(k, lo + i0 + )/ir. For QMC we obtain this function by the Maximum Entropy 
method from the dynamical Green's function G(k, r), while in our MFLEX approach we use the Pade approximation 
to continue the data from the line lo + iS slightly above the real frequency axis. Fig. |3| shows the data for U — At and 
T = 0.33i, where QMC is contrasted. This figure displays a very good overall agreement between both technique, 
thus strengthening our confidence in the MFLEX approximation. A similar comparison of the quasiparticle dispersion 
between MFLEX and QMC, but for the three-band Hubbard model was carried out in Ref. |3^. We now turn to the 
magnetic properties of the Hubbard model. Fig. ||(a) shows the static magnetic susceptibility Xzz(<\, w = 0) calculated 
within the MFLEX approximation [Eq. (]§)] compared with QMC data for U/t = 4 and different system sizes along 
the standard path in the Brillouin zone (0,0) — > (n, 0) — * (tt, tt) — > (0,0). We first consider this quantity for the 
comparison since it does not rely on an analytical continuation to the real frequency axis performed by the Maximum 
Entropy method and therefore the QMC result has very small errors. In Fig. |](a) we show the QMC results for 
U = 4t, x = 0.08, t' = — 0.2i, t" = 0.15i, T = 0.33t obtained on a 8 x 8 lattice. For the same parameter set and system 
size, the corresponding MFLEX results (squares) show a smaller susceptibility and a much smaller peak at (tt, tt). 
This means that the antiferromagnetic fluctuations are underestimated within the MFLEX approximation compared 
with those in the QMC simulations. The question arises whether this is due to the fact that we have neglected 
vertex corrections or, rather, to the approximation itself. Dahm et alEi determined the leading contributions to the 
vertex corrections (up to 0(U 2 )) and found that these tend to reduce the spin susceptibility. On the other hand, 
the comparison of our results with QMC apparently shows that the spin-spin correlation function calculated without 
vertex corrections are underestimated with respect to the (in principle exact) QMC results. Thus, either the quality 
of the perturbative results worsens when including vertex corrections or higher order contributions tOpthese diagrams 
become important at our intermediate values of U . In the case of a non self-consistent RPA calculationc3 it was shown 
that the introduction of an effective U* < U restored a good comparison with QMC calculation. This is because in 
the RPA approximation without self-consistent Green's function magnetic fluctuations are overestimated by the RPA 
denominator. In our calculations, x is renormalized self-consistently by the renormalization of G, which, in turn, 
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reduces \ especially near its maximum at (7T, 7r). However, this reduction overcompensates for the enhancement due 
to the RPA denominator and this is the reason why magnetic fluctuations are underestimated in this self-consistent 
calculation. Taking example from Ref. one could think to introduce an effective U* greater than U in order 
to compensate for the reduction of magnetic fluctuations. This procedure would also be a simplified version of the 
"pseudopotential approach" introduced by Bickers et al.c3 to include parquet diagrams in an effective way. However, 
we have verified that increasing U from At to 8t at a fixed high temperature T = 0.33t results in a slight decrease 
of Xzz(l)- The opposite occurs at lower temperatures like T = 0.05i. Thus the introduction of a temperature- 
independent U* (which is needed to extrapolate QMC results to lower temperatures) cannot improve the results for 
X- 

On the other hand, since magnetic fluctuations are very sensitive to the temperature, we introduce an effective 
temperature T* < T and compare our MFLEX results calculated with the temperature T* with QMC results calculated 
at the temperature T. The physical motivation for this ansatz is that, due the closeness of the system to the Mott- 
Hubbard transition at half-filling, antiferromagnetic fluctuations are very strong!!! These fluctuations, while fully 
captured by QMC, are underestimated by the MFLEX approximation which is not able to describe the metal-insulator 
transition appropriately. Since antiferromagnetic fluctuations are very sensitive to the temperature, the shortcoming 
of the MFLEX approach can be removed by a reduction of the temperature. Indeed, our comparison with QMC results 
for U — At is greatly improved if one takes a scale factor R such that T* = T/R with R = 2.5. The MFLEX results 
at this temperature T* can be seen in Fig. ||(a) (triangles) to compare quite well with QMC data at temperature T. 
This is especially true for the correlation length as indicated by the arrow. From the same figure, one can also infer 
the importance of using a method which allows to increase the system size. Here indeed, we also present the MFLEX 
results for a 64 x 64 system at T = 0.33i. For this increased system the static susceptibility shows significantly more 
pronounced magnetic fluctuations than the 8x8 results with the same T. It is also interesting that the 64 x 64 
MFLEX results resolve the small peak between (0,0) and (w, 0) observable in QMC although the magnetic response 
is still smaller than the one obtained with QMC at the same temperature. 

The same renormalized temperature T* can be used in order to achieve a good agreement with QMC of the 
imaginary part of the dynamic spin susceptibility Imx zz (Q,w). In Fig. ^(b) we show Imx^z(Q,w) for U = At 
obtained with MFLEX and QMC on the same small 8x8 systems. As for the static susceptibility, the MFLEX 
approximation yields a much smaller value for Imx 2Z (Q, oj) than the QMC result, whenever calculated with the same 
temperature T. However, if the MFLEX results are calculated at the reduced temperature T* introduced above 
the agreement is drastically improved. In particular, the slope of Imx zz (Q,w) at u> = (related to T\, see below) 
as well as the position u> s { of the maximum agree very well. Thus, renormalizing the temperature in the MFLEX 
calculations leads to a considerable improvement of the perturbative results as compared with QMC. We have verified 
that the effective temperature T* is related to the true temperature T by approximately the same scale factor R 
also at higher temperature. For T = t we obtain an optimized agreement with QMC results similar to Fig. |] with 
R « 2. This gives us confidence that the renormalization factor will be appropriate to describe qualitatively the spin 
correlation function also for smaller temperatures at least within an error of ~ 20% which we anyway allow for the fit 
to experiments. Finally, the temperature renormalization does not affect the agreement of the quasiparticle dispersion 
as shown in Fig. since the latter depends only weakly on the temperature. The importance of the temperature 
for the spin fluctuations becomes clear in Fig. B, where we show the (logarithm of) the inverse of the frequency ui s f 
where Inr\; Z z(Q, w) is maximum, as a function of U for different T. I/w s f is an indication of the strength of the 
antiferromagnetic fluctuations, since it diverges in the spin-density- wave state. For high temperatures, l/w s f is only 
weakly {/-dependent while for very low._temperatures uj s f vanishes exponentially with increasing U indicating the 
SDW instability (at T = 0) for U -> ocO Below, we will show that our results with U — At, (whereby the MFLEX 
temperature is set to the renormalized temperature T* = T/Rt), agree quite well with experimental results on T\ 
and reasonably well with T2G at low temperatures. 

Nevertheless, 4t is believed that the properties of the cuprate materials at and close to half filling are better described 
by a larger f//nL3. We will thus also show the results for U/t — 8. However, applying perturbation theory when the 
interaction is as large as the bandwidth (U = 8t « W) is questionable. Nevertheless, it is tempting again to compare 
the diagrammatic results with QMC and thus use the MFLEX calculation to extrapolate the QMC data to lower 
temperatures and larger system sizes. The static susceptibility for U = 8i on 8 X 8 clusters is presented in Fig. ||(a). 
For T — 0.33t, Xzz(q) obtained in the present approximation is rather flat and structureless and does not compare 
well with QMC results. We thus use again our strategy of renormalizing the temperature by a factor i?x, which must 
necessarily depend on U and is expected to be larger for increasing U. As one can see from Fig. ^|(a), we need a 
temperature renormalization factor i?x of about 5 in order to have a good agreement for the static spin susceptibility 
Xzz(q)- On the other hand, from Fig. ||(b) one can see that it is not possible to find an appropriate temperature 
renormalization which makes the imaginary part (i.e. the dynamical properties) Imx zz (Q, to agree with the QMC 
result. If one requires that only the slope at u) = coincides (which is the important quantity necessary to calculate 
Ti) we need Rt ~ 3. This value of Rt does not coincide with the one obtained for the static correlation function. This 
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clearly shows the difficulty of using this diagrammatic approach for such large interaction strength. For comparison 
with experiments in Sec. |V|, we will use an intermediate temperature renormalization factor J?t = 4. 

V. COMPARISON WITH EXPERIMENTS 

Most of the available experimental results on the magnetic properties of the high-T c materials are extracted from 
nuclear magnetic resonance (NMR) and inelastic neutron scattering (INS) studies. While INS measures directly the 
q and w-resolved Im% zz (q, u), NMR, as a local probe, determines weighted averages of the susceptibility Xzz(q_,u) 
for to — > over the whole Brillouin zone: Specifically, the spin-lattice relaxation time Ti probes the inverse of the 
slope of Imx zz (q, to) for to — > and the spin-echo decay time T-iq the inverse of the static susceptibility. NMR and 
INS experiments on LSCO and YBCO have revealed a lot of remarkable properties: (i) strong antiferromagnetic 
fluctuations persisting in the normal as well as in the superconducting phase up to the optimally doped regime, (ii) 
a suppression of Im% zz (q, to) at small to attributed to a spin gap opening at low temperatur.es in metallic YBCO, and 
(Hi) a sharp resonance peak at 41 meV and q = Q for optimally doped YBCO below T c .E§E-3 The spin gap manifests 
itself in INS measurements with a depression of the magnetic response at low energies and low temperatures. NMR 
measurements agree with this spin gap and show-a depression of T\T below T^ s which is about T c for the overdoped 
and larger than T c for the underdoped samples .E3 

To relate the spin-lattice relaxation time Ti to the spin susceptibility \zz i we adopt the approach by Shastry, Mila 
and Ricalj describing the hyperfine coupling of the Cu 2+ spins with the different nuclei in the unit cell, which leads 
to the expressions 

q 

where the form factor T c (q) results from the Fourier transform of the hyperfine interaction 

F c (q) = {Aa b + 2B[cosq x +cosq y ]} 2 . (12) 

Here, we consider the case where the applied static magnetic field is perpendicular to the Cu02 planes. A different 
geometry of the experiment would require a different form factor. 

The transverse relaxation rate T2G is related with the static spiiususceptibility through the Gaussian component of 
the spin echo, as pointed out independently by Thelen and PinescZl and TakigawacS: 



T -2^ 0-69 
2G 128U 2 



I]TF c 2 ff (q)xL(q,0) 



N 
q 



(13) 

with a different form factor F e ff(q) for T2G, obtained form Eq. ( p_2[ ) by replacing A a b with A c . 

The unknown hyperfine coupling constants A a b, A c , B are extracted from Knight shift experiments. Here, we adopt 
the values recently given in the analysis by Barzykin and PinesQ and set A a b = 0£AB-. A c = —AB and finally the 
energy scale B = 3.82 x 10 _7 eV. These are similar to values given by other authorsoL3iJ Note, that both relaxation 
times give complementary information: while T\ probes the slope of the imaginary part of Xzz(q, to) for lo — > 0, T%g 
depends on the static susceptibility x zz (q, to = 0). Since NMR probes the local environments of the spins, all momenta 
q contribute in principle to the relaxation rates as can be seen from Eqs. ( pl]Jl^ ). However, in the presence of a large 
antiferromagnetic correlation length, the q points close to the AF point (n, tt) will give the largest contribution to 
these expressions. 

In the following analysis, we choose U — 4.5i which turns out to give the best agreement with the experimental 



results regarding magnetic properties. Moreover, as discussed in Sec.tV, for intermediate U (U ^ 6t) only is it 
possible to have a good comparison with QMC results with a unique temperature-renormalization factor i?x- For 
larger U (U > 7t) this cannot be made unambiguously. The energy scale t is fixed by taking t = 250 meV, so that 
the bandwidth is ~ 2eV as observed in photoemission experiments. Furthermore, we take the same temperature 
renormalization factor Rt = 2.5 as found for U = At, since Rt does not change much from U — 4t-to U — 4.5i. The 
experimental resultsJbr Ti and T2G of the two YBa2Cu307_a samples are taken from Imai et al.E3 for YBa2Cu307 
and from TakigawaEj for YBa2Cu30@.63 (these data are collected in Ref. 0). 
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In Fig. 0(a) , we show the spin-lattice relaxation time T\ multiplied by T as obtained from the MFLEX calculations 
on a 64 x 64 lattice for x = 0.08 and x = 0.16 corresponding to the two samples YBa 2 Cu3 06.63 an d YBa2Cu3C>7, 
respectively. The figure shows that T\T is proportional to T for high temperatures, while for lower temperatures it 
tends to a constant. While the linear behavior at high temperatures agrees with experiments, for T smaller than 
a certain value, T\T should increase again due to the occurrence of a spin gap. In our calculation we are not able 
to see this gap behavior, possibly because we cannot reach the temperature where the gap sets in or because of 
the limitations of our approximation. However, a precursor of the spin gap is seen in the flattening of T\T at low 
temperatures. The linear behavior of T\T with temperature (indicating T\ = const.) aw.ay. from the spin gap regime 
is actually very well reproduced by MFLEX calculations in a wide range of parameters.OEil Aiming at carrying out 
a semi- quantitative comparison with experiment, it is thus natural to fit the two parameters of the linear behavior, 
namely, the T — » extrapolation and the slope of T\T versus T. The theoretical and experimental values obtained for 
these two parameters are listed in Tab. |[ Notice that we took into account the temperature renormalization factor 
i?T and modified the slope accordingly. In agreement with the experimental data, we find that the extrapolated T\T 
is only slightly larger for the YBa2Cu30y sample than for the underdoped YBa2Cu30g.63 sample. Moreover, the 
bigger slope present in the overdoped sample suggests that the two functions cross at some higher temperature T croS s 
(« 50K), as observed in experimental results. 

The 63 Cu spin-echo decay time T2G calculated according to Eq. ( |l3|) is shown in Fig. @(b). Again, the measured 
data show approximately a linear- T behavior in the range between ~100K and 300K, in agreement with our theoretical 
results. For a semi-quantitative comparison with experiments we again extract the slope and the T — > extrapolation 
and show the results in Tab. Note, that increasing the hole doping results in a shift of T2G to larger values, while 
the slope remains almost the same, in agreement with the experimental findings. 

While for T{T both the slope and the T — > 0-extrapolated values agree quite well with the experiments (within 
25%) in the case of T2G, the agreement is good only for the slope, while the T — > extrapolated values is too large, 
especially for the underdoped sample. In principle, we could try to adjust this extrapolated value by decreasing \t'\ 
but this would worsen the results for T\T . To show that a deviation of 25% is a good result, we consider the effect of 
a small change in t' . We thus include the data for U — 4.5<, t' — — 0.25t, t" — 0.15t in Tab. | showing that a change 
of 25% in t' results in more than 100% changes in the T — » extrapolated values of T\T and T2G- Notice that the 
slopes are not very much affected by such a change. By increasing U or decreasing \t'\ the values of T\T and T2G 
extrapolate to smaller and eventually to negative values. This signals that the system approaches a SDW instability. 
In Tab. § we also include results for U — 8t, but with t' = — 0.4, t" = 0.15 and a temperature renormalization factor 



of Rt = 4 as discussed in Sec. IV. The results are worse than the U — 4.5i ones, especially concerning the T — > 
extrapolation. Notice that the latter are quite sensitive to t' and could be improved by increasing \t'\. On the other 
hand, the slope is essentially independent on t' and cannot be improved in the same way. 

Notice that the too large value for the U = 4.5t, T — > 0-extrapolated Tig is due to the fact that this quantity is 
extremely sensitive to the value of U. For example, for U = 4t one would have obtained Tig ~ 200/is for T — > 0. A 
fine tuning of U could fix the extrapolated value of Tig more accurately, although this would put T\T off. 

Although the value U/t = 4.5 seems to give the best agreement with the magnetic properties at finite doping, the 
same value of U ft does not reproduce correctly the insulating behavior for the half filled model. Electron energy 
loss and optical experiments have revealed a charge transfer gap of ~ 1.7eV for YBCOlIj which would require rather 
large values for U/t > 8 for op: chosen value of t = 250meV. Moreover, INS experiments on the antiferromagnetic 
parent compound YBa2Cu30gL3 showed that this system is well described by a spin- 1/2 antiferromagnetic Heisenberg 
model with an exchange coupling of ~ 0.125eV. Using t = 250meV and J = 4t 2 /U one needs a value of U/t « 8. 
On the other hand, our diagrammatic approach cannot be well reliable for such a large value of U/t as discussed in 
Sec. IV, even when using a temperature renormalization factor extracted from the comparison with QMC data at high 
temperatures. For this reason, we cannot rule out that a more appropriate calculation could give a good comparison 
with experiments on T\T and Tig also for U/t ~ 8. Another possibility could be that the effective U/t at finite doping 
may be reduced with respect to the one at half filling due to the screening of the doped carriers. 

We now consider the correlation length of these systems. This can be inferred from the q-dependence of the 
static spin susceptibility x*z(<l,u = 0), plotted in Fig. |. for U = 4M,t' = -0.2, t" = 0.15, a; = 0.08, T = 150K. 
Since Xz Z (q, u> — 0) is strongly peaked at the antiferromagnetic wave vector Q = (jr, tt), the susceptibility appears 
to be commensurate. This is in contrast-with LSCO which clearly shows maxima at the incommensurate points 
Q' = (tt + 9o,tt) (and symmetric points) O'Ej Experimentally, it is not clear whether an incommensuraJ3ility is seen 
in the spin response of the YBCO materials. Early INS experimentsEJ and some more recent onesE3 suggest a 
commensurate structure, while other authors report experimental data that are better fitted bv-a superposition of 
Lonentzian curves peaked at the four equivalent incommensurate points (tt, tt ± q ) and (tt ± q , 7r) E3 or at (tt ±q ,n± 
q c )J£3 where the incommensurability q Q is material and doping dependent. To extract the correlation length £ and the 
incommensurability q Q of the calculated spin-spin correlation function we model Xzz (q, u — 0) by the superposition 
of four Lorentzian curves with width l/£ peaked at the points (tt, tt ± q Q ) and (tt ± q a , tt). We carry out the fit along 
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the line {%, q) with < q < n where the model function discussed above reads 

r 

Xzz(q,u = o) = 



(l/C) 2 + {q-Q + q ) 2 
r 2r 

(14) 



(i/O 2 + {q-Q- qo) 2 (i/O 2 + (« - Q) 2 + <z 2 

with Q — tt. The results for the correlation length are shown in Fig. ^|. As observed in experiments, the correlation 
length £ is of the order of 1 — 2 lattice spacingsoE3o and temperature independent for T < 250K. It decreases 
for increasing doping levels. Although the single maximum in the curve of Fig. ^| suggests a commensurate struc- 
ture, it is better fitted with an incommensurability q a « 0.5. This agrees with recent INS experiments —where an 
incommensurability of q a w 0.4 was suggested to better reproduce the spin susceptibility in YBa2Cu306.6-£3 

If we define £ to be just the half width half maximum (HWHM), the correlation length is even smaller (£/a » 1 
for x — 0.08) but still temperature independent. Our finding of a relatively small correlation length £/a ~ lg — 2 is 
in contrast with the phenomenological NAFL treatment of the NMR relaxation times by Barzykin and PinesJH where 
rather large correlation lengths of about 7a for YBa2Cu3 06.63 and 2a for YBa2Cu30y were necessary for satisfying 
fits of the experiments. That these correlation lengths are too large in. comparison with experiments was already 
pointed out in a later critical reexamination by Zha, Barzykin and Pines.Ej A relatively small correlation length is also 
the reason why the temperature dependence of £ and T2G & r e different. 

Finally, in Fig. [l^ we plot as a function of T the magnetic energy scale u> s { defined to be the energy where lm.Xzz (q, 
takes its maximum. The calculated values are [between 20 and 40 meV, in good agreement with the energy of the 
magnetic resonance peak found in experiments.cil A comparison between Fig. ^ and Fig. [h] suggests that T\T and 
u>si show the same linear temperature dependence, in agreement with the NAFL theory. For T — » 0, w s f also tends to 
a constant which decreases with decreasing \t'\ or increasing U, i. e., approaching the SDW instability. 



VI. SUMMARY AND CONCLUSIONS 



In summary, we have studied the electronic and magnetic properties of an underdoped and an overdoped 
YBa2Cu307_,5 sample with S = 0.37 and 5 = 0, respectively. We started from the two-dimensional Hubbard model 
including longer-ranged hopping processes to describe the correlation effects in these materials. Since it is essential 
to reach both low temperatures and a fine spectral resolution for a qualitative comparison with experiments, we em- 
ployed the MFLEX approximation, which amounts to neglecting particle-particle fluctuations and vertex corrections 
in the FLEX approximation. We checked the quality of this approach by comparing with Quantum-Monte Carlo 
results at higher temperatures. We found that for not too large Hubbard interactions U (U < W/2) the agreement 
between the MFLEX and QMC results can be considerably improved by introducing a [/-dependent renormalized 
temperature T* = T/Rf. We then use this temperature renormalization factor i?x in order to extrapolate the QMC 
results towards the temperature necessary for an analysis of the experiments. One should be aware of the fact that 
this approach is uncontrolled and may be ineffective, since the temperature renormalization factor Rt found at the 
temperatures accessible by QMC (T/t ~ 0.33) may not hold for lower temperatures and thus the extrapolation may 
fail. On the other hand, as we already mentioned, the temperatures accessible to QMC are far too high and thus an 
extrapolation scheme to lower temperatures is mandatory. 

In the search for appropriate parameters of the Hubbard model to describe qualitatively and semi-quantitatively the 
magnetic properties of YBa2Cu30y_5, we find a fair agreement with experimental results on T\ and T2G, (within an 
average error of less than 20%) for the parameter set U = 4.5t, t' = — 0.2<, t" = 0.15t and t = 250meV. For this value 
oiU/t, we need a temperature renormalization R w 2.5, as inferred from the comparison with QMC. Our calculations 
describe in a qualitative way also the shape of the Fermi-surface, and the flat quasiparticle energy dispersion near 
(tt, 0) which approaches the Fermi surface at optimal dopingcJ. Moreover, we have a reasonable description of (i) the 
slope and T — > 0-extrapolated values of T\T, where T\ is the spin lattice relaxation time, (ii) the slope of the spin-echo 
decay time Tig vs - temperature, (iii) the correlation length £ w 1-2 lattice spacings, (iv) the q-dependence of the 
spin response, which appears commensurate due to a single maximum at Q = (tt, tt), but which is better fitted by 
a superposition of four Lorentzian curves peaked at incommensurate peaks, and (v) the typical size of the magnetic 
excitation energy scale of 20-40 meV. Finally, our calculation is not inconsistent with a parameter set with stronger 
coupling, like U/t = 8, as seems to be the case for the cuprate materials, provided one increases the value of the 
next-nearest-neighbor hopping. However, our procedure is less reliable in this parameter regime. 

Our results thus indicate the importance of introducing finite values for longer-range hoppings t' jt and t"/t for the 
sake of a qualitative description of magnetic properties. This is in agreement with numerical calculations on the i tt J 
model which show the relevance of t'/t and t"/t for an appropriate description of single-particle properties. E3cilEJ. 
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Notwithstanding the importance of longer-range hoppings appearent from these works, it is difficult to establish 
whether other parameters, like possibly an interplane coupling or longer-range interactions may be physically more 
important and thus describe the experiments more appropriately even without introducing t' and t". 
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Slope 
[fl] 


T -» 
[sK] 


Slope 
[MsK- 1 ] 


T -> 

H 


YBa 2 Cu 3 7 
MFLEX x = 0.16 


5.0- 10 -4 
4.0 ■ 10" 4 


7.5 • 10" 2 
7.9 ■ 10" 2 


0.18 
0.26 


84 
153 


YBa 2 Cu 3 6 . 63 
MFLEX x = 0.08 


4.0 • 10 -4 
3.2 ■ 10" 4 


8.5 • 10" 2 
8.4 • 10" 2 


0.14 
0.19 


23 
124 


MFLEX x = 0.08 

U = 4.5i, t' = -0.25t, t" = 0.15t 


1.8 ■ 10" 4 


2.2 • 10 _i 


0.11 


198 


MFLEX x = 0.08 

U = 8t, t' = -0.40t, t" = 0.15* 


2.8 • 10~ 4 


1.9 • 10 -2 


0.24 


102 



TABLE I. Slope and extrapolated value for T greater than the spin gap obtained by fitting a straight line to the measuredLJLj 
(labelled with YBa 2 Cu 3 7 -,5) and calculated (MFLEX) data for T\T and T 2G with U/t = 4.5, t'/t = -0.2,t"/t = 0.15. 
Comparison with the results for U/t = 4.5, t'/t = -0.25, t"/t = 0.15 and U/t = 8, t'/t = -0.4, t"/t = 0.15 is shown. 
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FIG. 2. Quasiparticle dispersions for the Fermi-surfaces 
shown in Fig. [[]. We only show the most interesting flat part 
around (n, 0). 
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FIG. 3. Spectral functions A(k,uj) 

for U = At,t' = -0.2t,t" = 0.15*, a; = 0.08, T = 0.33t as 
obtained by QMC (top) and MFLEX (bottom) on 8 x 8 sys- 
tems. Dark (bright) areas correspond to a large (low) spectral 
weight. Note that the energy dispersion depends weakly on 
the temperature. 
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FIG. 4. (a) Static spin susceptibility x^A"\^ u = 0) 
along the standard path in the Brillouin zone as obtained 
by MFLEX and QMC calculations. (b) Corresponding 
dynamical spin susceptibility ImXzz(Q,<^>) at the antifer- 
romagnetic momentum Q = (tv,tv). All data are for 
U = 4t,t' = -0.2t,t" = 0.15t and a 8 x 8 lattice except 
for the MFLEX results in (a), additionally labeled by 64 2 , 
which are for a 64 x 64 system. 
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FIG. 5. The inverse of the magnetic energy 

scale u! s i as a function of U for various temperatures 
(x = 0.08, t' = -0.2, t" = 0.15). A SDW instability mani- 
fests itself by uj Bi — » 0. 
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FIG. 6. Same as Fig. |, but for U = 8t and x = 0.10. (a) 
Static spin susceptibility (b) Corresponding dynamical spin 
susceptibility 
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FIG. 7. (a) TiT and (b) T 2G NMR relaxation times 
calculated in the MFLEX approximation for U/t = 4.5, 
t' /t = —0.2, and t"/t = 0.15, and two different doping lev- 
els, representing YBa2Cu3 06.63 (a; = 0.08) and YBa2Cu3C>7 
(x = 0.16). 




FIG. 8. Static susceptibility in the Brillouin zone for 
U = 4.5t, x = 0.08 and T = 150K.Although this function is 
strongly peaked at Q = (n, n) indicating a commensurate spin 
response, it is better fitted by four incommensurate peaks. 
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FIG. 9. Correlation length obtained by fitting the static 
susceptibility along the line from (n, 0) to (n, n) with a super- 
position of 4 incommensurate Lorentzian-shaped peaks sym- 
metrically located around Q. 
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FIG. 10. The magnetic energy scale lo b { defined to be the 
frequency where the dynamical susceptibility at Q is maximal. 
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